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Abstract. We study the effect of radiation pressure on the dust in the inner rim of transitional disks with large 
inner holes. In particular, we evaluate whether radiation pressure can be responsible for keeping the inner holes 
dust-free, while allowing gas accretion to proceed. This has been proposed in a paper by Chiang and Murray-Clay 
(2007, Nature Physics 3, p. 604) who explain the formation of these holes as an inside-out evacuation due to X- 
ray-triggered accretion of the innermost layer of the disk rim outside of the hole. We show that radiation pressure 
is clearly incapable of stopping dust from flowing into the hole because of dust pile-up and optical depth effects, 
and also because of viscous mixing. Other mechanisms need to be found to explain the persistence of the opacity 
hole in the presence of accretion, and we speculate on possible solutions. 
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1. Introduction 

In recent years our understanding of the structure of the 
dusty disks around young stars has increased considerably 
with the enormous data volumes of such disks from space- 
based infrared telescopes (ISO, Spitzer Space Telescope), 
ground-based infrared interferometry, and millimeter wave 
observations. One interesting issue that emerged is that a 
certain fraction of these protoplanetary disks appear to 
have large inner holes. To be more precise, these objects 
have inner disks that are nearly dust-free out to several 
AU, much farther out than the dust sublimation radius. 
Such disks are frequently termed "transitional disks" , ex- 
pressing the idea that these disks are in the process of 
becoming debris disks. 

Several ideas for the origin of these holes are being 
discussed. One idea is that photoevaporation might have 



opacity holes without the usual dust content. Chiang & 



drilled a hole in the inner disk ( 


Hollenbach & Gorti 


2005 


Alexander et al.||2006 


Clarke 2007 Ercolano et al. 


2008) 



Another explanation is the presence of planetary systems 
inside the "gap" (e.g. Alexander & Armitage 2009). 



Surprisingly, several transitional disks still show signs 



of accretion (e.g. Fang et al. 2009 Sicilia-Aguilar et al. 



Murray-Clay (2007 henceforth CM07) propose a mecha- 
nism by which the ionizing effect of X-rays from the star 
causes an inside-out evacuation of the disk. In this sce- 
nario, radiation pressure on dust grains plays a crucial role 
in holding back the dust while allowing the gas to flow in. 
CM07 make a detailed analysis of the conditions needed 
to activate the magneto-rotational instability (MRI) nec- 
essary for accretion. In their model, the disk is largely 
viscously inactive because of insufficient ionization. Only 
directly at the inner rim of the disk. X-rays at energies of 
about 3 keV from the central star lead to sufficient ioniza- 
tion in a thin layer with column density ~ 5 x lO^'^cm"^. 
As this layer viscously spreads inward and is being ac- 
creted by the star, new layers get exposed to the ionizing 
radiation. In this way, the disks gets evacuated from the 
inside out, leading to a steadily growing inner hole. CM07 
find accretion rates that range from 10"^"'^ to 10~^ Mq /jr. 
If the dust was dragged along with the gas, and if the ac- 
cretion onto the star occurs in a continuous way through 
a disk, then there will be no opacity hole. CM07 propose 
that radiation pressure on dust grains is the mechanism 
to keep the inner disk free of dust. While CM07 consider 
the possibility that some dust leaks into the inner hole, 
there is no mechanism to eventually remove it again from 
the hole. 



2010 1, indicating that some gas is streaming through these 
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Here we revisit this issue of holding back the dust by 
radiation pressure. Our main concern is the optical depth 
effect, which limits the total amount of dust that can be 
affected directly by radiation from the central star. As 
the accretion layer proceeds into the disk, dust will accu- 
mulate and eventually collectively overcome the radiation 
pressure. 



2. Model 

Let us consider a radial column of gas and dust in the in- 
ner disk rim, exposed to radiation from the star. Radiation 
pressure acts only on dust grains in the inner few optical 
mean free paths (r* < few), with the radiation force de- 
creasing with increasing optical depth. Ignoring other ef- 
fects, the innermost grains would therefore quickly catch 
up with particles at higher optical depth, creating an in- 
finitely thin sheet of dust. Particle diffusion (either ther- 
mal or turbulent) will act to keep the thickness of this 
sheet finite. Furthermore, it will also transport the effect of 
radiation pressure to the dust deeper in the disk, thereby 
allowing a larger amount of dust to be affected. The radi- 
ation pressure will thus act as an outward moving "snow 
plow", piling up more and more dust. 

As the dust particles are pushed outward through the 
gas, they will feel a friction force. The total friction force 
on the pile will increase as the pile builds up. On the 
other hand, the radiation pressure is fixed. All particles 
in the pile will have to share this given amount of radia- 
tion pressure. The dust pile will therefore gradually slow 
down. Even a very small inward pointing accretion flow 
will eventually stop the snow plow and then drag the dust 
pile inward with the gas. 



2.1. Analytical estimates 

2.1.1. When does a pile-up take place at all 

The radiation force of the unattenuated stellar luminos- 
ity Li, acting on a single grain with radius a in the rim 
surface located at a distance ro from the star is given by 
-Frad = Li.Qpj-Tra'^ / {4:TrrQc) where Qp^ is the momentum 
absorption efficiency of the grain and c is the speed of 
light. For the analytical and numerical estimates below, 
we take Qpi- = 2, to simulate the most favorable condition 
for holding grains back with radiation pressure. At the 
same time, the gas is streaming slowly toward the star, 
with the accretion velocity Wacc- In the Epstein regime 
and for subsonic dust drift, we can write th e friction force 

c;t;th/07ra^/3 



per grain as F^: 



le.g. 



Schaaf 1963) 



where p is the gas density and Vth = \/8kT/ {-Kjimp) is the 
thermal velocity of gas particles with mass fxrup at a tem- 
perature T. rup is the proton mass, and we take fj, = 2.3 
for molecular gas. 

The first question we need to answer is whether there 
will be any pile-up at all. This is the case if the radiation 
force onto a single grain exceeds the friction force onto 



that same grain. The criterion for this is 
3 L-i,Qp,- 



7 := 



IGtT r^CTaccWthP 



> 1 



(1) 



If 7 < 1, the dust will accrete with the gas onto the star 
without delay; i.e., the space velocity of the dust grains 
will never point away from the star. If 7 > 1, on the other 
hand, the accretion of dust will be delayed until a large 
enough pile has accumulated. 

2.1.2. Speed of the snow plow relative to the gas 

For the moment we will assume that the gas in the outer 
disk is not moving at all, and that the radiation pres- 
sure is simply driving the snow plow a distance s into the 
disk, starting at a position tq. The speed of the column 
piled up by radiation pressure can be calculated by equat- 
ing the radiation force onto a column with 1 cm^ inner 
surface /rad = / (AnrQc) , with the inward-directed col- 
lective friction force /fric of all the grains in the pile while 
the pile is moving relative to the gas with a velocity Vdr- 
The latter force is /mc = o-piicWdrWthp/(^a), where cTpiie is 
the column mass density of dust in the pile and ^ is the 
specific density of the grains which we take to be 2g/cm^. 
Solving for the drift velocity we find 



47rr^C(TpiieWthP 



(2) 



2.1.3. Time-dependent solution and total mass 
supported by radiation pressure 

For simplicity we assume that the gas density p and the 
dust-to-gas ration /dg are constant prior to buildup of the 
pile. We define a variable s that measures the distance 
from the rim into the disk, so that we can write the posi- 
tion of the dust pile as rpi\e(t) = tq -\- s(t). Then the dust 
mass in the pile can be expressed easily as f7piie(s) = pfdgS. 
The drift velocity at location s into the disk is then 



Vdr{s) 



K 



(3) 



'i-Krlcp^ fdgVt\,s ' s 

where K :— L*^a/(47rrQcp^/dgWth)- As we can see, the 
velocity of the pile becomes lower and lower as it collects 
more dust to be pushed forward. We can integrate p]) t o 
find the distance traveled by the pile to be s{t) = \f2Kt. 

Now we introduce a constant inward gas motion with 
velocity Wacc- The position of the pile now becomes 
''piic(0 = ^0 + ^(t) — t'acci, while s(t) remains unchanged. 
The snowplow comes to a halt at tstop = -^/(^i^acc) ^-^d 
rstop = ro -I- K/v^cc - Wacc^stop- After that, continued mass 
loading will turn around the motion of the pile and move 
it, embedded into the accretion flow, into the inner hole 
and toward the star. The total column at the turn-around 
point is given by 

^acc 



f^pilc('5stop) 



(4) 
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Table 1. Some numerical estimates for observed transitional disks 



Star 


ivh 


I/* 


M 


a 




7 


Cpilc 


rstop/ro 


^stop 








[IQ-iOAfo/yr] 




[AU] 




[lO-^g/cm^] 




[yr] 


DM Tau 




0.34*^) 


10.0(1) 


0.035(1) 


3(1) 


0.69 


1.4 


1.00008 


0.6 


TW Hya 


0.6(1) 


0.23'^) 


4.0(1) 


0.01(1) 


4(1) 


1.19 


2.4 


1.0001 


4.3 


GM Aur 


1.2(1) 


1.33'^) 


50.0(1) 


0.007(1) 


25(1) 


0.97 


1.9 


1.00004 


15.3 


HD100546 


2.5(*) 


36'*) 


io.o(*) 


0.035(^) 


10(4) 


57.57 


115.1 


1.03464 


602.4 



References: (1) CM07; (2) Ricci et al. (20101; (3) Calvet et al. (20021; (4) Dominik et al. (20031; (5) copied from DM Tau 



2.1.4. Width of the peak from diffusion 

We can also estimate the width Ai? of the dust pile as it 
is produced by the interplay between accretion and vis- 
cosity. The diffusion time across Ai? is given by tdiff = 
AR^ /D, and the accretion time across the same distance 
is iacc = AR/v^cc- Equating these timescales, and using 
D = v/Sc where D is the diffusion coefficient, v is the 
viscosity (molecular or turbulent) and Sc is the Schmidt 
number, we find Ai? = v / {Scvs.cc)- Please see below for 
more details on the Schmidt number we use for our calcu- 
lations. 



2.2. Numerical values 

Table [1] shows numerical values for the quantities derived 
above. To derive these values we have relied on stan- 
dard accretion theory (e.g. Hartmann 2009[ ) and used the 
following relations: Wacc = 3i^/(2r), S = M /{2ttRvs,cc), 



and p — I]/(v27r x hp) with the pressure scale height 



hp '■= Cs/^Ik- 



The table shows 7, apHc, rgtop, and igtop for the three 
stars in CM07 and also for HD 100546, a Herbig star that 
also shows a large gap even though is does have some hot 



dust close to the star (Benisty et al. 2010). We have in- 



cluded this star because of its higher luminosity, which is 
helpful for illustrating the pile- up effect. The value of 7 
shows that the radiation force on individual grains does 
not greatly exceed the friction force created by the accre- 
tion flow - something one would have expected to achieve 
significant pile-up. Consequently, the total column of dust 
piled up before the turn-around point is only of the order 
of 1 X 10~^ g cm~^, and the location where this happens 
is no more than a few times 10~^ AU, very close indeed to 
the original starting point of the rim. The pile-up happen 
on a time scale between 1 and 40 years. So even starting 
with a clean inner disk, dust will start to flow into the gap 
after that time. 

The higher luminosity of the Herbig star HD 100546 
does have a slightly stronger effect, but it will be not any 
more noticeable than in the case of T Tauri stars. 

3. Numerical model 

3.1. Model description 

The calculations so far are based on a number of sim- 
plifying assumptions to make the problem analytically 



tractable. We also constructed a numerical model to study 
the radial dust motion in this situation. Instead of solving 
the complete coupled set of equations involving both the 
dust and the gas motion, our goal is less ambitious. We 
want to see whether the radiation pressure can counter 
the dust pile-up enough to keep the inner hole dust-free. 
Therefore, we study the reduced problem of radial 1-D 
dust motion with radiation pressure in a disk with a given 
inward gas motion. 

The model is built in the following way. We start by 
setting up the gas density Pg(r) at the disk midplane 
as a function of the radial coordinate r. The gas den- 
sity is computed semi-selfconsistently by specifying the 
accretion rate Al in the disk, the turbulence parame- 
ter for accretion aa and the flaring angle if for stellar 
irradiation. Using the flaring angle, the midplane disk 



temperature is estimated as (Chiang & Goldreich 1997): 
Tmid{r) = (/7^/'*(r*/r)^/^T^,, where i?* is the stellar radius 
and is the stellar effective temperature. We assume 
that no viscous heating takes place (valid for low accre- 
tion rates and large inner holes) and that the full 7ri?J 
stellar surface is "visible" at any point on the disk surface. 
From the resulting temperature we can now compute the 
isothermal sound speed Cs — y/kT/pnip, where k is the 
Boltzmann constant and fi ~ 2.3 is the mean molecular 
weight of the gas and rup is the proton mass. The accretion 
viscosity coefficient is now given as Va = aaC^/^K, where 
flK = y/GAU/r^ with the stellar mass. Finally, with 
AI = ^TTYigVa (valid for r ^ R^,) the gas surface density 
can be computed for given value of the accretion rate 
AI. For more details on standard accretion disk theory. 



we refer to the book by Hartmann ( 2009 ) . We do included 
neither the effect of the dust motion nor of changing dust- 
to-gas ratios, in our estimation of the gas disk structure. 

We then insert the dust into this model. We assume 
that the dust grains are small enough that they do not 
settle much, so that the increase in midplane dust den- 
sity due to settling is negligible. We focus solely on the 
radial motion of the dust. At each time step we start with 
a very simple radiative transfer calculation: We follow the 
radiation of the star as it moves along the disk equatorial 
plane outward, i.e. we integrate the equation^ "^^^ ^^^^^ = 
—pd{r,t)KdL{r), where pd{r,t) is the dust density as a 
function of spatial position r and time t. is the Planck- 
mean dust opacity, computed with the stellar temperature 
as reference temperature. For simplicity we assume that 
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Fig. 1. Results of the numerical model for the source 
HD100546 for the case Sc=l, i.e. A = i^t and a* = 0.035. 
Top panel: Surface density of dust as a function of radius 
(horizontal axis) and time (vertical axis). The color is the 
dust density in normalized units. The contour levels of the 
diagram are steps of 10°'^, i.e. five contours are one factor 
of 10. The level in the bottom-right corner is the initial 
value of the dust surface density, 2.3 x lO^'^ g/cm^. The 
bottom panel is the same, but now plotted as different 
curves for different times, with the vertical axis being the 
dust surface density. The model actually calculates the 
density of the dust pd, but for the plotting we converted 
it to surface density = y/2TTHpPd. In this model the 
curves that are more to the top are the later ones. The 
time interval between two curves is 10 years. 



Kd is independent of r, i.e. we neglect the changes in stellar 
spectrum as it is extincted in the disk. In fact, we do not 
solve for the frequency-dependent luminosity L^(r,t), but 
instead for the frequency- integrated one L{r,t). This ex- 
erts a radiation pressure force on the dust grains, in units 
of dyne/gram-of-dust, of fr{r,t) = KdL{r,t)/{Ancr'^). The 
dust has friction with the gas and we compute the dust 
velocity by assuming that radiative and friction forces are 
in equilibrium, i.e. that the dust drifts relative to the ac- 
creting gas at its equilibrium drift velocity. 

In addition to this systematic dust motion, there is 
also random motion due to the gas turbulence. While tur- 
bulence may be extremely weak in a dead zone, we shall 
see that even the slightest turbulence will profoundly af- 
fect the structure of the dust pile-up region. The turbulent 
"diffusion" of dust is due to the same turbulence that is 
believed to be responsible for the accretion process (the 




0.0010 



0.0001 



Fig. 2. As Fig. [l] but now for Sc=103, i.e. A = 0.001 ft. 
Time differences between curves is 300 years this time. 

turbulent "viscosity" i/a, see above). The turbulent diffu- 
sion constant Dt is therefore of the same order of magni- 
tude as i>a - Exactly what the ratio Da/va is is not perfectly 



understood yet (Johansen & Klahr 2005). If viscosity is 



due to turbulent gas motions, this ratio should be not 
far from unity. However, it may also be that large scale 
magnetic fields could remove angular momentum from the 



disk without inducing much turbulent mixing ( Combet & 
Ferreira 2008). Here we parametrize it as the Schmidt 



number: Dt = vt/Sc. We treat Sc as a free tuning param- 
eter. 

The dust motion now follows from the continuity equa- 
tion: 



dpd{r,t) 1 d{r'^pd{r,t)vd{r,t)) 



dt 



]_d_ 

j,2 Qj. 



dr 



d ( Pd{r,t) 



dr \ pg{r) 



(5) 







which we solve using implicit integration ("backward 
Euler"). 

As an initial condition for the dust we set pd to 0.01 pg 
for r > To and = for r < tq, where tq is the initial 
dust inner rim radius, given in table [ij We take the dust 
grains to be a — 0.1 pvn in radius, consisting of material 
with a density of ^ = 3 g cm~'^. For simplicity we take 
a simple geometric opacity oi k — 3/(^a), which is twice 
the geometric cross section. As boundary condition at the 
inner edge of the computational domain we take a simple 
outflow condition, while at the outer edge of the domain 
we set Pd = O.Olpg as a boundary condition. 
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Fig. 3. As Fig. [ij but now for the source DM Tau. Time 
interval between curves is 200 years. The contour levels of 
the top panel are steps of lO"'^, i.e. five contours are one 
factor of 10. The level in the bottom-right corner is the 
initial value of the dust surface density, 1 x 10"^ g/cm^. 



3.2. Model results 

We performed the numerical equivalent models for all 
sources listed in table [T] One of the most striking results, 
found in all models, is that if we take Sc=l, i.e. Dt = vt, 
then any outward movement of the dust induced by the 
radiative pressure is completely overwhelmed by the enor- 
mously strong turbulent diffusive mixing. So even for cases 
where the radiative force would easily push out the grains 
if no turbulent mixing would occur (Sc— >■ oo), the grains 
are washed inward as a result of strong mixing. An exam- 
ple of this is shown in Fig. [l] In this example the time 
scale for mixing dust toward the inner edge of the com- 
putational domain at 9.5 AU is about 10 years only. The 
overall solution quickly becomes dominated by the outer 
boundary condition (best seen in the bottom panel), ren- 
dering it meaningless at that point. 

Why is particle mixing so dominant, almost indepen- 
dent of at? This can be seen in the following way: If, we 
take strong/weak turbulence, e.g. a — Q.l / a = 10^^, 
then, for a given accretion rate M, the gas surface den- 
sity becomes low/high (by virtue of M = STrEz^t). So for 
a given radiation pressure force on the grains, the drift 
speed is high/low, but also the turbulent mixing as well 
as the inward gas velocity become high/low. Working out 
the numbers for the accretion rates given in table [T] in all 
cases the turbulent diffusion wins over radiation pressure. 



The conclusion is that with Sc ~ 1, the radiation pressure 
is entirely unable to prevent the dust from drifting inward 
toward the star. 

However, as discussed earlier, if angular momentum 
transport is not primarily caused by turbulence, much 
larger Schmidt numbers are allowed. We choose Sc=10^, 
but keep the at the same as before. Then we find that 
for some cases, notably HD 100546, the radiation pressure 
is strong enough to push the dust outward against both 
the accretion flow of the gas and the action of turbulent 
diffusion. For other cases the gas drag of the inward gas 
motion is still too strong, such as for DM Tau. The results 
for HD 100546 and DM Tau are shown in Figs.([2) [3|. 

For HD 100546 one can see that the radiation pressure 
pushes the grains that are in plain sight of the star out- 
ward. But the absorption of the radiation (the same pro- 
cess that transfers radiative momentum to the grains) also 
shields the dust grains at larger optical depth from stellar 
light. These dust grains therefore still flow inward with the 
gas. What happens is, as one can see in the flgure, that 
the dust piles up in a rim of high dust density. According 
to our model, if one takes an even higher Schmidt num- 
ber this dust pile-up rim becomes geometrically thinner. 
However, the total mass in this rim (and thus its optical 
depth) is only dependent on the amount of dust that has 
been swept up. The radiation force acts only on the inner- 
most layer of this dust pile-up rim, i.e. the layer of radial 
optical depth unity. But the diffusion, even if ever so weak, 
will distributed the momentum that is received by these 
inner dust grains to the entire mass of dust in the pile-up 
rim. One also clearly sees that the dust pile-up rim first 
moves outward, but then slows down as more and more of 
the radiation momentum has to be divided over an ever 
increasing mass of dust in the pile-up rim. This dilution of 
the radiation force means that the outward motion of the 
rim stalls, and eventually the dust gets dragged inward 
with the gas flow. 

For DM Tau we see that, even with a high Schmidt 
number, a pile-up never really happens, because the gas 
drag by the inward moving gas is simply too strong. One 
sees a nearly instant inward motion taking place, as ex- 
pected based on the 7 value in Table [T] being lower than 
unity. 

4. Discussion and conclusions 

In summary, wc flnd strong arguments that radiation pres- 
sure alone is wholly insufficient to keep the inner hole of 
a transitional disk dust-free while low-level accretion is 
continuing. Both dust pile-up and viscous mixing are eas- 
ily able to overcome the radiation pressure quickly and 
refill the inner hole with dust on a level similar to nor- 
mal dust-to-gas ratios. These processes pose a problem 
for the inside-out evacuation scenario for transitional disk 
in its original form. While we do think the MRI activa- 



tion discussed by Chiang & Murray-Clay (2007) may be 



responsible for launching a low-level accretion from the 
inner rim, mechanisms other than radiative forces have to 
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be invoked to keep dust from flowing into the inner holes, 
or at least from contributing to the opacity in these holes. 
An additional force to hold back dust that has been 



discussed in the literature is photophoresis (Krauss & 



Wurm 2005), which is a strong force caused by the tem- 
perature difference between the illuminated and dark sides 
of a dust grain. This force may have an effect for larger 
grains (cm-sizes), but it is irrelevant for the sub-micron- 
sized grains considered here (see figure 1 from Krauss & 



Wurm| ( [2005| )). 

One possible solution is that the dust has grown to 
such large sizes in these regions that the dust continuum 

A 



2005) 



optical depth drops below unity(Tanaka et al. 
way to achieve this could be efficient dust growth in the 
pile-up produced by radiation pressure In fact, the jump in 
gas surface density will foster faster and more efficient dust 
settling, further enhancing the dust density. Greatly in- 
creased coagulation rates might be the result. Large grains 
produced in this way could even enter the gap without be- 
ing visible, and without being an obstacle to the ionizing 
radiation needed to drive the inside-out evacuation of the 
disk. However, a much more detailed model of this sce- 
nario is required before any conclusions can be drawn. 

A second possibility is that the inner hole has been 
made by a planet orbiting close to the disk rim, and that a 
pressure bump created by the planet-disk interaction traps 



Rice et al. 


2006 


). However, 


Ward 



(2009) has questioned the effectiveness of this mechanism. 

A third possibility would be that torques exerted by 
one or more planets located in the inner disk hole in the 
system transport matter so quickly to the star, that the 
average optical depth in the inner disk remains small even 



though dust is contained in the accreting matter (Zhu 



et ah 2011 Perez-Becker & Chiang 2011). 



Our model contains a number of simplifications, for 
example the use of geometrical cross sections for the dust 
grains. We also did not account for the optical depth of 
the gas inward of the dust rim. The gas is expected to 
have much lower optical depth than the dust, but at high 
enough accretion rates even the gas can become optically 
thick to stellar radiation (Muzerolle et al. 2004). This 
effect is rather weak at the low accretion rates we are 
concerned with (according to Muzerolle et al. it plays vir- 
tually no role at M ~ lO~^M0/yr), and therefore it is not 
problematic that we ignore it. As with other assumptions 
in this paper, we have always made choices that would 
make it easier to hold back the grains by radiation pres- 
sure. The conclusion that dust grains cannot be held back 
by radiation pressure is therefore a firm one. 
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